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We investigate the maximal violation of Bell inequalities for two d-dimensional systems by using 
the method of Bell operator. The maximal violation corresponds to the maximal eigenvalue of the 
Bell operator matrix. The eigenvectors corresponding to these eigenvalues are described by asym- 
metric entangled states. We estimate the maximum value of the eigenvalue for large dimension. A 
family of elegant entangled states |^)app that violate Bell inequality more strongly than the maxi- 
mally entangled state but are somewhat close to these eigenvectors is presented. These approximate 
'\^'' states can potentially be useful for quantum cryptography as well as many other important fields 

' of quantum information. 

o : 

, PACS numbers: 03.65.Ud, 03.67.-a,42.50.-p 

a : 

^ ^ I The famous Clauser-Horne-Shimony-Holt (CHSH) inequality 0,01 for two entangled spin-1/2 particles has always 
provided an excellent test-bed for experimental verification of quantum mechanics against the predictions of local 
^ ■ realism |^. It is well-known that all pure entangled states of 2 dimension violate the CHSH inequality: the maximum 
quantum violation of 2-^/2 being often called the Tsirelson's bound 0. In 2000, contrary to previous study, Kasz- 
^ likowski et al showed numerically, based on linear optimization techniques, that the violations of local realism increase 
t — , with dimensions for two maximally entangled d-dimensional systems (qudits) (3 < d < 9) 1^. A year later, Durt et al 
■ extended the analysis to d = 16 under special experimental settings|^. In the same year, an analytical proof was con- 
. structed for two maximally entangled 3-dimensional systems (qutrits) 0. In 2002, two research teams independently 
developed Bell inequalities for high-dimensional systems: the first one is a Clauser-Horne type (probability) inequality 
, for two qutrits ^ ; and the second one is a CHSH type (correlation) inequality to two arbitrary c?-dimensional systems 
' 0, now known as the Collins-Gisin-Linden-Masser-Popoescu (CGLMP) inequalities. The tightness of the CGLMP 
. inequality was demonstrated in Ref. 1^. The maximally entangled state of two-qudit reads |^')mos — X]j=o \^^) 
C^, where |j) is the orthonormal base in each subsystem. Collins et al restricted their investigation to |^')nios, because 
on one hand |^')mos is a simple state, on the other hand there is a very natural viewpiont that maximally entangled 
^ , states would maximally violate the Bell inequality, just as the CHSH inequality has worked for two-qubit. The results 
^ ' of was numerically consistent with Ref. For |^')mcs, the Tsirelson's bound or Bell expression Id is given by 



X 



^ 01, 1 

/,(!*)_) = 4d^(l-^)(^ 



,^0 - ^ 2d3sin2[^(fc + i)/d] 

^' " 2d3sin2[7r(-fc-l + i)/d]^' 

showing that the maximum quantum violation increases with the dimension d (here I — [d/2— 1] represents the integer 
part of [d/2 — 1). ). When d — > 00, it is interesting to note that 

2 °° 1 1 
lim /d(|*),„cs) ^ — Y.' 



7r2 ^^(fc + 1/4)2 (fc + 3/4)2J 
~ 2.96981, (2) 

a number related to the Catalan's constant. 

Contrary to prevalent belief, Acin et al studied the quantum nonlocality of two-qutrit as well as two d-dimensional 
systems up to d = 8 and discovered another unexpected result: there existed non-maximally entangled states that lead 
to greater violation of the CGLMP inequalities compared with maximally entangled states 0|. This surprising result 
still lacks of a good intuitive explanation. From the TABLE I of Ref. 11], one observes that the maximal violation 
increases with dimension d, and it reaches 3.1013 for d = 8. This gives rise to a natural question:" What are the 
maximal violations of the CGLMP inequalities if one increases the dimension d, especially when d goes to infinity?" 
The purpose of this Letter is to investigate the maximal violations of the CGLMP inequalities for higher dimensional 
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systems. We also present a family of elegant entangled states |^')app whose corresponding maximal violations are 
approximate to the real ones. 

Before comput ing the quantum violation, let us first estimate the upper bound under quantum mechanics. Sig- 
nificantly, Ref. [13 has recast the CGLMP inequality into a form that is very similar to the CHSH inequality: 

/d=Qll+Ql2-Q21+Q22<2, (3) 

where Qij are the correlation functions defined by probabilities in the following way 

^ d-l d-l 
m— n— 

where S — (d — l)/2 is the spin of the particle for the d-dimensional system, /'■' (m, n) — S — M {e{i — j){m + n), d); 
e(x ) = 1 and —1 for x > and a; < respectively; M{x,d) — {x modd) and < M{x,d) < d—1. Refs. and 
[l3| have proved that Id < '2 for hidden variable theory. From Eq. one notes that the extreme values of Qij are 
±1 for both local realistic description and quantum mechanics, therefore it is impossible that the maximum quantum 
violation of 1^ is larger than 4. Furthermore cannot reach 4, because Qij are constrained to each other, if three of 
them are set to be 1, the fourth must also be 1. Consequently, one can conclude easily from above analysis that the 
maximal quantum violation of is a number between 2 and 4. 

To generalize Acin's approach, we first note that it has been shown that unbiased multiport beam splitter |l2l| can 
be used to test violation of local realism of two maximally entangled qudits. Unbiased d-port beam splitter is a device 
with the following property: if a photon enters any of the d single input ports, its chances of exit are equally split 
among the d output ports. In fact one can always build the device with the distinguishing trait that the elements of 
its unitary transition matrix T are solely powers of the root of unity 7 = exp(i27r/d), namely Tki = "^7'^'- In front 
of i-th input port of the device a phase shifter is placed to change the phase of the incoming photon by 4>{i). These d 
phase shifts, denoted for convenience as a "vector" of phase shifts (j) — {4>{0) , 4>{^) , ■■■,4>{d — 1)), are macroscopic local 
parameters that can be changed by the observer. Therefore, unbiased d-port beam splitter together with the d phase 
shifters perform the unitary transformation U{4>) with the entries Uki = Tki exp{i(j){l)) . Devices (Bell multiports) 
endowed with such a rnatrix were proposed, and readers who are interested in it can refer to Refs. IHIT^. The 
approach developed in 11] is related to Bell operator. An arbitrary entangled state of two-qudit reads 

d-l 

The quantum prediction of the joint probability P{Aa ~ k, ~ I) when Aa and _Bf, are measured in the initial state 
1$) is given by 

P{Aa=k,Bb = l) = Tr[(C/((^„)t®C/((^fc)t)rife®ri, (f/(0a)®C/(<P6)) |*)(*|] 

d-l 

— \^ n ,n* pil<l>aU)+VbU') + ^{jk-j'l)~4'a{m)-ipb{m')-^{m.k-m'l)] (f.\ 

j^j' ,m,m'— 

where 11^ — \k){k\, IT; = are the projectors for systems A and B, respectively. Substituting Eq. ^ into the 
CGLMP inequality, one gets the Bell expression for the state 1$): 
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y^^f,il'l>lU)-'l>l{m)+VlU')-Vlim')] - )(g'¥fc(.7-m) _ ^-i 2- (fc+1) (j' -m') ^ ^ 

t 91. 

gi[02O)-<!i2(m)+V2(i')-¥'2('n')] _ ) (g» ^ fcQ-m) _ ^-^jliL (fe+1) (j' -m') -j| _ 

can be expressed as 



/d(|$))=Tr(B|$)($|) = (<i>|B|$), (8) 

where B is the so-called Bell operator. Startmg with the CGLMP inequality and choosing suitable the experimental 
settings |^[Tlt 

optimal for maximal violations, we can derive the element of the Bell operator matrix as 



d-l 



1=0 

k=a ^ ^ 

« 01, 

eT0-™)+55(/-™') ^(1 f^)(e-'if ('=^+1)0'-™') _ e*if /^O'-™)) + 

(i-™)-!?!/-™') I^)(e*if '^O-™) _ g-^if (^+1)0'-™'))}^ (10) 

fc=o ^ "'^ 

where B, in general, is a x (P matrix. Ref. |l j| found that the maximal eigenvalue of matrix B is nothing but 
the highest quantum prediction of the Bell expression Idi\^) cig) and the corresponding eigenvector |<&)oig is the state 
that maximally violates the Bell inequality. Thus, the problem of computing the maximal violation of /^(|$)) reduces 
to the determination of the maximal eigenvalue of matrix B. Due to X)f=o^ e^^ip^i)^ = dSpq, where 6pq = 1 when 
p = q modulo d and otherwise, the matrix B can be further simplified, i.e., it can be decomposed into the sum of d 
decoupled operators that act individually within the subspaces spanned by the vectors {|00), |11), \{d— 1))}, 
{|01), |12), ...,\{d- l)0)},...,{|0(d- 1)), |10),..., |(d - l){d - 2))}, respectively. For the first subspace spanned by the 
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vectors {|00), |11), \{d — l){d — 1))}, and with the constraint j — m = j' 



m', one can reduce Bmm'.jj' to 
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= - sm 
a 
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fc=0 
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{j - to) 
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d^l 

1 



) sin 



d — 1 COs[^{j — TO 



r(l - (^mj), 



(11) 



redl 



where B^^^ is the matrix element located at the TO-th row and the j-th column of the reduced Bell operator B 
Let us denote = B^.m+r, since -Bmj depends only on the difference (j — m). Therefore, we arrive at the following 
d X d real symmetric matrix for the reduced Bell operator: 



B 



redl 
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Bl 


B2 
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Bd-i 




Si 


Bo 


Bl 


B2 




Bd-3 


Bd-2 




B2 


Bl 


Bo 








Bd-3 






B2 








B2 






Bd-3 








Bo 


Bl 


B2 




Bd~2 


Bd-3 




B2 


Bl 


Bo 


Bl 




{Bd-1 


Bd-2 


Bd~3 




B2 


Bl 


Bo 


) 



(12) 

One may read from Ea. (|ll|l that B^ are all non-negative numbers, Bo = 0, and interestingly liuid^oc Bd-i = 4/7r, 
lim(i^oo Bd-j / Bd-i = for j << d. The Bell operator B'^^'^^ acts on a matrix column (oq, ai, • • • , a^-i)"^, which 
corresponds to the state l^*) J2'j=o the first subspace. Actually, an arbitrary state |$) = J2'j^^=o'^n'\^j') 

can always be transformed into its Schmidt decomposition form |^) = X]j=o through local unitary transforma- 

tions, thus it is sufficient to study the maximal violation problem in the first subspace. For a general state |^), one 
can obtain the Bell expression as 



d-l d-1 



rn— j— 



d-l . d-l-r 
r—1 ^ m— 



(13) 



For instance, for the maximally entangled state |4')rnes with aj — l/Vd, the summation (X]r?x j=o -^'"j)/'^ ^ 
(Sr=i 2('^ — r)Br)/d recovers the results of Eq. 

As the analysis shown in the maximal violation of the CGLMP inequality with the experimental settings 
@ corresponds to the maximal eigenvalue of B'^^''-^^ and indeed its corresponding eigenvector is a nonmaximally 
entangled state of two-qudit. For instance, for d = 3, one has Bq = 0,Bi = 2^/3/3, B2 = 2, the eigenvector 
YY3^(|00) + v^-v^ |ll) + |22)) corresponds to a maximal violation /d=3(|^')cig) = 1 + ^11/3 - 2.9149, 

2.8729; for d = 
1 



which is larger than /d^ad^E' 
Bl 



4, one has Bq = 0,Bi = 2^4-2x72/3,52 = 2\/2/3, 



2V4 + 2V2/3,the eigenvector |*; 



(i=4 
oig 



= (|00) + a\n) + a\22) + |33)), with a = (V2 + V2 + 



8 - 3V2 + 4^2 - \/2- V4 + 2V2) /{V2+ V4^-2V5) ~ 0.73937, corresponds to a maximal violation /<i=4(|*)oig) 



+ + I V 8 - 3V2 + 4\/2- V2 ~ 2.9727, which is larger than /d=4(|*) 



2.8962. 



By diagonalizing exactly the matrix B'^^^.'^^ we have obtained the real maximal violations Id{\'^)cig) for two entangled 
qudits. The highest dimension that we have calculated is c? = 8000. In Fig. 1, one may observe that /rf(l^')cig) 
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FIG. 1: Variations of J^d'!')) with increasing dimension d {2 < d < 8000). 



increases with dimension d. Based on the data of Id{\'^)cig) from d = 2 to d = 8000, one has an empirical formula 
fitting /d(l^')cig) numerically to the dimension d: 

C''''(l*>cig) - 3.9132- 1.2891x-"-2226 (14) 

from which one can see that /™"*^^(|*)cig) ^ 3.9132 is a coarse-grained limit of the maximal violation for the CGLMP 
inequality when d tends to infinity. 

Analysis of the eigenvectors |^')eig shows that these eigenvectors numerically satisfy some general properties: for 
instance, |\E')cig — J2jZo with maximal eigenvalue has the following symmetric properties for the coefficients: 

flj = ad-i-j', and oq : fli : a2 : 03 : • • • ~ 1 : -i= : -i= : -j= : ■ ■ ■ for large d. Thus, we may approximate a family of 
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elegant entangled states 



I*) 



,app 

'j 



1 1 



app — 



V77 ^(j + i){d-j) 



1 



(15) 



{j + l){d-j) 



whose corresponding Bell expressions /(i(|5')app) are closed to the actual ones /d(|^)cig)- For example, for d ~ 8000, 
the error rate between /d(l^')cig) and /d(|^)app) is only about 0.745%. We have also listed /ci(l^')app) and /d(l^')mcs) 
in Table I and drawn the corresponding curves in Fig. 1. One may observe that for d — 50000, /(i(|5')mcs) — 2.96981 
has almost reached the limit as shown in Eq. |21); and for d — 600000, /c;(|^')app) has exceeded 3.80. 

In summary, we have investigated the maximal violation of the CGLMP inequalities for two entangled qudits. The 
maximal violation occurs at the non-maximally entangled state, which is the eigenvector of Bell operator with the 
maximal eigenvalue. The maximal violations increase with growing dimension d, coming close to 4 when d approaches 
to infinity. Experimental test has been performed to verify the CGLMP inequalities for the first few high-dimensional 
systems, and indeed there exist non-maximally entangled states that violate these inequalities more strongly than the 
maximal entangled ones |il4 | . Bell inequalities are applicable to quantum cryptography and quantum communication 
complexity |l5| . previous researches are mostly based on maximally entangled states. For stronger violation, it may 
be useful to generate the approximate state in this paper. Nevertheless, it may be significant and interesting to apply 
the elegant symmetric entangled states |^')app to quantum cryptography as well as other important fields of quantum 
information. 
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d 


2 


3 


4 


5 


6 


7 


8 


9 


10 


20 


30 


40 


50 


60 


/d(l*)cig) 


2.82843 


2.9149 


2.9727 


3.0157 


3.0497 


3.0777 


3.1013 


3.1217 


3.1396 


3.2492 


3.3068 


3.3449 


3.3728 


3.3946 


Jd(l*>app) 


2.82843 


2.90909 


2.96466 


3.00689 


3.04078 


3.06895 


3.09296 


3.10639 


3.13219 


3.24551 


3.30496 


3.34389 


3.37222 


3.39417 


/d(|*)mes) 


2.82843 


2.87293 


2.89624 


2.91054 


2.92020 


2.92716 


2.93241 


2.93651 


2.93980 


2.95472 


2.95974 


2.96225 


2.96376 


2.96477 


d 


70 


80 


90 


100 


150 


200 


250 


300 


350 


400 


450 


500 


550 


600 


/d(|*)cig) 


3.4124 


3.4273 


3.4400 


3.4511 


3.4914 


3.5178 


3.5370 


3.552 


3.5641 


3.5743 


3.583 


3.5906 


3.5973 


3.6033 


/d(l*>app) 


3.41192 


3.4267 


3.4393 


3.45022 


3.48933 


3.51446 


3.53256 


3.54651 


3.55775 


3.5671 


3.57505 


3.58194 


3.588 


3.59339 




2. !)():) i!) 


2.9GG():i 


2. !)()() 15 


2.9()()7S 


2.9G779 


2.9()8:!0 


2.9G«G() 


2.9()880 


2.9G895 


2.9()90G 


2.9G9il 


2.9()92i 


2.9G92G 


2.9()9:!1 


d 


650 


700 


750 


800 


850 


900 


950 


1000 


1500 


2000 


2500 


3000 


3500 


4000 


/d(|*)cig) 


3.6087 


3.6136 


3.6181 


3.6222 


3.6260 


3.6296 


3.6329 


3.6360 


3.6594 


3.6747 


3.6859 


3.6946 


3.7017 


3.7077 


^ci(|*>app) 


3.59824 


3.60263 


3.60664 


3.61032 


3.61372 


3.61687 


3.61981 


3.62256 


3.64299 


3.65623 


3.66584 


3.67330 


3.67936 


3.68442 


/d(|*)meB) 


2.96935 


2.96938 


2.96941 


2.96944 


2.96946 


2.96948 


2.96950 


2.96951 


2.96961 


2.96966 


2.96969 


2.96971 


2.96973 


2.96974 


d 


5000 


6000 


7000 


8000 


50000 


70000 


80000 


90000 


100000 


200000 


300000 


400000 


500000 


600000 


/d(l*)cig) 


3.7174 


3.7250 


3.7311 


3.7362 






















^d(l*>app) 


3.69253 


3.69884 


3.70398 


3.70829 


3.75659 


3.76372 


3.76644 


3.76878 


3.77083 


3.78345 


3.79019 


3.79472 


3.79810 


3.80080 


/d(|*>mes) 


2.96975 


2.96976 


2.96977 


2.96978 


2.96981 


2.96981 


2.96981 


2.96981 


2.96981 


2.96981 


2.96981 


2.96981 


2.96981 


2.96981 



TABLE I: [Optional Table (for referee)] Bell expressions /d(|^')eig), /d(|*)app) and /d(|*)mes)- 



